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We introduce templates for exponential asymptotic expansions that, in contrast 

to matched asymptotic approaches, enable the simultaneous satisfaction of both 

boundary values in classes of linear and nonlinear equations that are singularly 

perturbed with an asymptotic parameter e — > 0"'' and have a single boundary layer 

j^^l at one end of the interval. For linear equations, the template is a transseries that 

takes the form of a sliding ladder of exponential scales. For nonlinear equations, the 

transseries template is a two-dimensional array of exponential scales that tilts and 

^^ . realigns asymptotic balances as the interval is traversed. An exponential asymptotic 

^Sj ' approach also reveals how boundary value problems force the surprising presence 

of transseries in the linear case and negative powers of e terms in the series beyond 

all orders in the nonlinear case. We also demonstrate how these transseries can 

be resummed to generate multiple-scales-type approximations that can generate 

uniformly better approximations to the exact solution out to larger values of the 

r^ . perturbation parameter. Finally we show for a specific example how a reordering of 

^ ' the terms in the exponential asymptotics can lead to an acceleration of the accuracy 

of a truncated expansion. 
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'nI I In this paper we examine the role of exponential asymptotics in boundary value 

problems (BVPs) for classes of second-order ordinary differential equations with a 
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'^ ' small positive asymptotic parameter e multiplying the highest derivative. In par- 
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ticular we consider both linear and nonlinear problems with boundary layers at 
one end of the finite domain. At first sight, this appears to be a well-trodden 
classical problem with little more to be revealed. Systems that cannot be solved 
exactly can be attacked using a vast toolbox of asymptotic approaches, most no- 
^^ . tably matched asymptotic expansions (MAE), see for example Nayfeh (1973) and 

5h ' references therein. 

^ Despite, or perhaps because of, the success of MAE approaches, BVPs have re- 

ceived less attention from asymptotic analysts than initial value problems. Results 
on existence and uniqueness of solutions of second-order nonlinear two-point BVPs 
predated MAE (Wasow 1956 and references therein). Formal theorems on the ex- 
istence and character of asymptotic solutions of BVPs have been proved (Wasow 
1956; Howes 1978). More recent asymptotic BVP work (e.g., O'Malley & Ward 
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1997; Ou & Wong 2003, 2004; Ward 2005) has considered shock fronts and multiple 
spike solutions. However, we shall show that, even for such a well-studied problem, 
the deployment of an exponential-asymptotic approach throws up subtle, surprising 
and generic asymptotic issues left unanswered by previous approaches. 

Exponential asymptotics (Berry 1989; Segur et al 1991) seeks to obtain a for- 
mal expansion of a function that incorporates all relevant exponential scales. These 
scales are manifested typically by the exponential prefactors of (often) divergent 
infinite series. Techniques exist for managing the divergence and analytical con- 
tinuation of these series (Berry & Howls 1991; Olde Daalhuis 1993), which lead 
to exponentially improved numerical approximations that are uniformly valid in 
wider ranges of the underlying parameters (Olde Daalhuis 1998, Olde Daalhuis & 
Olver 1998). This approach can provide better error bounds (Boyd 1994) and al- 
gebraic methods to resolve complex geometric structures (Howls 1997; Delabaere 
& Howls 2002). In this paper we shall not be concerned with the implementation 
of the hyperasymptotics. Rather, we seek to identify the templates, or form for 
the expansions that are required as input to these hyperasymptotic approaches. 
The templates are the atomic building blocks which can be manipulated, using 
resummation techniques, to produce more sophisticated, 'molecular asymptotic' 
approximations. 

Exponential asymptotic approaches to solutions of linear and nonlinear ordinary 
differential equations have been well explored (Chapman et al 1998; Olde Daalhuis 
1998, 2004, 2005a&; Howls & Olde Daalhuis 2003; Olafsdottir et al 2005), but have 
concentrated on initial value problems. O'Malley and Ward (1997) applied expo- 
nential asymptotics to linear convection-diffusion PDE BVPs. Lee & Ward (1995) 
considered exponentially sensitive ill-conditioning in nonlinear ODE BVPs. Ward 
(2005) also provided an asymptotic description of multiple spike solutions and the 
associated interior boundary layers in Carrier-Pearson-type systems. However this 
paper is the first time that a detailed analysis of the satisfaction of singularly per- 
turbed BVPs at each and every exponential order has taken place. The approach 
reveals the following surprising results. 

First, as is well known, MAE usually violate one of the boundary conditions. 
This error may be negligible at small values of the asymptotic parameter e, but 
as e grows the range of validity (both numerical and analytical) of the expansion 
may decrease. Here we show how, for a similar effort, an exponential asymptotic 
approach satisfies both boundary conditions and hence obtains an expansion that 
is valid in a wider range with at least a comparable accuracy. 

Second, an exponential approach also reveals a generic, yet previously unknown, 
intricate and interlocking ladder (linear) or array (nonlinear) structure within the 
expansions that is essential to the satisfaction of boundary conditions. 

Third, for initial value problems arising from linear ordinary differential systems, 
usually only a finite number of asymptotic series play a role, each being a leading 
order (in the exponential-asymptotic sense) expansion of one of the solutions. The 
analysis here demonstrates explicitly how boundary values force the presence of 
"transseries" in linear systems here being infinite sums of exponentially-prefactored 
asymptotic series in the small parameter e. Transseries occur naturally in nonlin- 
ear differential systems (Olde Daalhuis 2005 o6; Costin 1998; Costin & Costin 2001; 
Chapman et al 2007) due to the mixing of exponential scales by nonlinear terms, 
but their role in linear systems has not been highlighted explicitly before. They are 
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significant, because the resummation of transseries can generate non-local informa- 
tion about the singularity structure of the problem and can be used to extend the 
range of validity of solutions. 

Fourth, for nonlinear BVPs, exponential asymptotics uncovers the presence of, 
at first sight, paradoxical negative powers of the small parameter e in expansions. 
Fortunately these are contained within series that are premultiplied by decaying 
exponentials in 1/e and so behave regularly as e — >■ O"*". 

Fifth, there is a choice in the order in which terms are resummed and this may 
be exploited, where possible, to allow for the derivation of numerically more efficient 
approximations to the solutions of nonlinear BVPs. 

At the outset, we stress that we are not proposing that an exponential- asymptotic 
approach should replace that of MAE (or any other) in boundary layer calculations. 
MAE remains a powerful and (where appropriately applied) numerically accurate 
approach. We aim to highlight the underlying asymptotic structure of BVPs to 
facilitate extensions to a matched asymptotic approach when needed. 

The outline of the paper is as follows. In section 2, we study a pedgagocial 
linear BVP to highlight the issues from an exponential asymptotic viewpoint. In 
section 3 we use general arguments to introduce the presence of transseries and the 
associated sliding ladder of exponential scales, demonstrating in section 4 how this 
will be a generic phenomenon in linear BVPs. In section 5 we introduce nonlinear 
BVPs, again with a pedagogical example, and the transseries array structure that 
allows boundary value satisfaction at both ends of the interval. The coefficients in 
the leading exponential orders are derived in sections 6 and 7 with a demonstration 
of the intricate realignment of balancing of terms, together with a demonstration 
of how it is possible to resum them to generate multiple-scales solutions. Moving to 
higher order terms in section 8 we uncover a contradiction that forces the presence 
of negative powers of the small parameter in the transseries. In section 9 we discuss 
how the reordering of terms within transseries can accelerate numerical agreement. 
We end with a discussion and suggestions for future work. 

2. A linear example of boundary failure 

We take a pedagogical approach and first illustrate the issues we wish to address 
through a specific example. The general case will be explained later. We consider 

eu"{x) + {2x + l)u\x)+2u{x)={), < a; < 1, (2.1) 

M(0) = a, u(l)=/3, 0<e«l, (2.2) 

where prime denotes differentiation with respect to a;. Conventionally the argument 
for a matched asymptotic approach runs along the following lines. Substitution of 
the standard ansatz 

oo 

u{x)r^^ar{xy (2.3) 

into (|2.ip generates first order linear differential equations for ar{x). Hence ar{x) 
cannot in general satisfy the two boundary conditions simultaneously. A boundary 
layer exists and there is an apparent need for a second expansion that obeys a 
rescaled equation, valid in the boundary layer and which satisfies the inner boundary 
condition. The inner expansion is then matched to the outer expansion involving 
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ar{x) that satisfies the other boundary condition. The inner and outer expansions 
are matched to determine the remaining unknown constants. 

Rescahng (j2.ip with a change of variables x = eX we may derive a composite 
expansion (at leading order) 

Mmatchcd(a;) = 3/3/(1 + 2x) + (a - 3/3)e-"/^ (2.4) 

Clearly we have Winatchcd(O) = a. However u(l) — j3 + O (e^^/') and so fails to 
satisfy the boundary condition at a; = 1, by at least the size of the inner expansion. 
For small values of e, this may well be a numerically negligible error, but this error 
will grow with e, so reducing the range of validity of the solution. 

In this example, we learn more by following the approach of Latta (Nayfeh 1973, 
pp. 145-154). We try the simple exponential asymptotics approach of a WKB ansatz 

oo oo 

u{x) - Y^ a^(a;)e'' + e-^^^^Z" ^ 6,.(2;)e^ (2.5) 

r=0 r=0 

We shall require F{0) — Q so that terms in the two series can balance at 0{e^) when 
x = 0. However, for general x ^ 0, 1, substitution of (12.51) into (12.11) and balancing 
at each order 0{e^) and O (e"^^^-'/'^), we obtain: 



0{e°): (2a; + l)ao(x) + 2ao(x) = 0, (2.6) 

0{e^): {2x + l)a'^{x) + 2ai{x) = -a'^{x), (2.7) 

O (^e-^(")/^e-i) : F'{xf - {2x + l)F'{x) =0. (2.8) 

The last equation coupled with F{Q) =0 gives 

F{x)=x^+x. (2.9) 

F{x) is positive on < a; < 1, as it should be, by assumption, so that (at least 
from an exponential-asymptotic viewpoint) no turning point occurs on that interval 
that can give rise to an interior boundary layer. The remaining equations for a^ 
and hr can be solved order by order. The simplified equations for br{x) are given 
below. 

The boundary conditions at x = are 

a^(0) + &r(0) == 5r^a, (2.10) 

where 6ij is the Kronecker delta. These can be satisfied by design. However the 
boundary condition at x = 1 can still only be satisfied up to exponential accuracy 
in e: 

a,(l)+o(e-^(i'/^) -^,0/3. (2.11) 

Ignoring the exponential error at the right-hand boundary we generate the relations: 

a,(x) = (a;_i(l)-a;_i(x))/(2x + l), ao(x) = 3/3/(2x + 1), (2.12) 

b'^{x)^h';_^{x)/{2x + l), hr{Q) = -ar{Q), bo{x)^a^3l3. (2.13) 
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The latter equations show that the hr are here actuahy constants. Hence the leading 
order solution is 



UWKB 



{x) = (ip/{2x + 1) + (a - 3/3)e-("'+")/') (1 + O (e)) . (2.14) 

Note that for a; > 0, the exponentially small term could be absorbed into the 0(e) 
term, but is kept separate to delineate the existence of the boundary layer. 

The matched and Latta-WKB approaches provide asymptotic "solutions" that 
satisfy the boundary condition near to the boundary layer at x = 0. They both fail 
to satisfy the boundary condition at a; = 1. From an exponential asymptotic (and 
for x{x + l)/e = 0(1), a numerical) viewpoint, this is unsatisfactory. Comparison 
of (|2.4p and (j2.14p shows that, for a similar amount of effort, we have Umatchod(l) = 
(^ + (a - 3/3)e"iA) (i + Q {e)) , uwkb(1) - (/? + (a - 3/3)e-2A) [l + O (e)) . The 
order of the exponential error in the WKB approach at the right-hand boundary 
is here less than in the matched case. This suggests that an extension of the WKB 
approach may be worth pursuing. We develop this approach in the next section. 

3. A transseries approach 

The exponential prefactors of the a^-series and fc^-series in (j2.5p balance at x — 
by design. At x = 1, the a^ satisfy the boundary data, but the br do not. 

Now suppose there is an additional series present, satisfying the differential 
equation (|2.1I) . but prefactored by e"^^^-'/'^. At x = 1, this can be used to cancel 
off the contribution from the fo^-series. However, while addressing the satisfaction 
of the boundary condition at x = 1 , the presence of this new series means that the 
boundary condition at x = is now violated exponentially. 

This new series can still satisfy the boundary conditions at x = if there is 
another, additional, series prefactored by e~'^^'^^'>+^W)/<^^ Yov then, at x = with 
-F(O) = 0, this second new series balances the first new series at 0(e^^^^-'/*). 

In turn, this second series can satisfy the boundary condition at x = 1 if there 
is another, third additional series prefactored by e~^^'^'' '^. In turn this can satisfy 
the X = condition if there is another, fourth additional series prefactored by 
g-(F(x)+2F(i))/e g^^ gp pj-^ g^^^ gp forth. (Problems with a single boundary layer 

near x = 1 would require F{l)—Q in (|2.9p and a reversal of roles of F{Q) and 
F{\).) In summary, the boundary values generate a ladder of series with exponential 
prefactors defining the runs (see figure [1]). 

The middle column in figure [T] can be visualised as a sliding scale of prefactors 
that can be used to line up different exponential scales. At x = the scale lines up 
with both with the 0-prefactor and also the scales at pF{\), for integer p > (left 
hand panel). This ensures that the boundary conditions at x = are satisfied. As 
X travels from to 1 this scale slides down (middle panel) so as to line up with the 
exponential prefactors pi^(l), p > 1 at x = 1 (right-hand panel). The terms that are 
algebraic in e in the re-balanced transseries can then be paired up so as to satisfy 
the boundary conditions, as we show below. With this ladder, the zero-prefactored 
series exactly satisfies the boundary condition at x = 1, with no exponential error. 

The template for the full asymptotic expansion is thus: 

oo oo oo oo 

ut.a„s(x; e) ~ ^ e-^^^^y^ ^ a^J>\x)e^ + ^ e-(^(-Hpi=^(i))/. ^ ^ir>) ^^y ^ (3.1) 

p— r— p— r— 
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Figure 1. Ladder structure of exponentials within the transseries ((3TTJ- The on the right 
of each panel denotes the series H2.3|) . The p_F(l),p = 1, 2, 3 . . . on the right-hand side of 
each panel denotes the series prefactored by exp(— p_F(l)/e). For general x the transseries 
exponents F(x) -|-pF(l), p = 0, 1, 2, 3 . . . on the left translate up and down so that at each 
boundary there is a balance available at each exponential order. 



where a}, (x) = ar{x). Such a sum of exponentially-prefactored series is called 
a transseries (Olde Daalhuis 2005a, 6; Costin 2001; Chapman et al 2007). The p- 
sums are transseries with prefactors that only depend on e, not x. When approached 
directly from the asymptotics, the presence of transseries in linear equations is, at 
first sight, surprising. 

We shall show below that this is the general form of the exponential-asymptotic 
template for the class of problems under consideration. More general templates 
might be derived for higher order equations. The template for BVPs with internal 
boundary layers must be modified to take account of the associated connection 
problems. 

Substituting (j3.ll) into (j2.ip and hereafter using -^^ to denote cc-derivatives, we 
arrive at the recurrence relations 



4rf(x) = [ai^\jl)^3b(r''>(l)-ai''\jx))/{2x + l), (3.2) 

(3.3) 



6(^(x) = bl^\^^{x)/{2x + l). 



^(-1) 



with 6}. '{x) = 0. The BVP then becomes 

aiP\0) + biP\0) = SroSpoa, 4'')(1) + biP-^\l) = Sro6po(3. 
Simple calculations give the leading orders of each series as 



(3.4) 



,(0)/ 



Ap) 



ip), 



al'>{x)=il3/{2x + l), a)^'\x) = -ZP{a-?.l3)/{2x + l), blf'{x) ^ S^ia - 3/3). 

(3.5) 
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Hence we have 

oo "^ 

+ (a - 3/3) J2 3Pe^(^(")+P^(i»/'= i (1 + 0{e)). (3.6) 

p=0 J 

The p-sums are convergent for |3e~^^"'^^/'^| < 1. Summing these is equivalent to 
changing the order of the p and r sums in transseries template p.ip . (Olde Daalhuis 
2005a, b; Costin 2001; Chapman et al 2007). The resuh is 

ut (x-e)-('^-/^^i^2:^^UI^^^^^le-™^Vl + 0(e)) 
utranslx, ^J " (^ 3^ + 1 \ 1 - 3e-^(i)/^ J 1 1 " Se'^^)/^ j j ^^ ^'' ■ 

(3.7) 
The approximation (J3.7I) satisfies both boundary conditions. In this sense (at 
least) approximation p.7p is better than either of the WKB or MAE approaches. 
If terms in e~-^^^^^'^ are neglected in p.7p . we recover the WKB result (|2.14p . 

Higher order approximations can in principle be derived by similar resumma- 
tions of the transseries at each order of 0{e'^). The results satisfy the boundary 
conditions at each order. In Appendix A we show how the p-series may be re- 
summed by a multiple scales approach. 

In figure [2] we compare the leading order results of the exponential asymptotic, 
Latta-WKB and matched asymptotic approaches with the exact solution for a = 1 , 
/3 = and the comparatively large small parameter e = 1. Clearly the exponential 
asymptotic approach is the only one to satisfy both boundary conditions. For the 
large value of e = 1, the other asymptotic approaches (especially MAE) agree 
noticeably less well with the exact result away from the boundary layer near a; = 1. 



»■" ^matched ,.--'' 

"~ ^ ^matched .. ' ' 




■-.MWKB 




Figure 2. A comparison of the exact and various asymptotic solutions of 1)2. 1|1 for a = 1, 
l3 — 0, e = 1 over the range < a; < 1 (left). The bold line is the exact solution, the 
thin solid line is the exponential asymptotic approximation, the dashed line is the WKB 
approximation, the dot-dashed line is the leading order MAE. On the right the integrated 
leading order error 1)3. 8p over the range < a: < 1 as a function of < e < 1 is plotted. 



In the right-hand graph in figure [2] we compare the integral of the relative errors 
over the entire range < 2; < 1: 

Aapprox(e) = / |wcxact(a;;e) -Uapprox(a;;e)|dx, (3.8) 

JQ 
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for values of e between and 1. Clearly the cumulated error in the exponential 
asymptotic approach is smaller over the entire range of x. 

4. Comparison with the exact solution and the general case 

A similar analysis can be carried out for more general solutions of second-order 
linear two-point boundary value problems. Consider the system: 

eu"{x)+c{x)u'{x)+d{x)u{x) ^0, < x < 1, (4.1) 

u{0) = a, u(l) = /?, <£<<!, (4.2) 

where c{x),d{x) > for < x < 1. Let the solutions of the equation be 

ui(x; e) = Tt\e), u^ix- e) ^ c-^(-)/%("\e), (4.3) 

where the T (e) and F{x) are given by 

N-l i-x 

T^'^\e)=J2f^i^e^ + R,,^{<:,e), F{x) = c{rj)d^. (4.4) 

For a choice of ui and U2 satisfying (|4.1I) . the self-consistent solution of the BVP is 

_ [0^2(1) - /3m2(0)] ui{x) + [/3mi(0) - auijl)] U2ix) 

"^^''^" [u,{o)Mi)-Mi)Mo)] ' ^ ' 

provided the denominator does not vanish. This can be written as 



1 _ e-^(i)/^Tf '(e)Tf )(e)/ri('\e)Tf )(e) j rf ^(e 
P _ ae-^(i)/^Tf ' (e)/rf ^ (e) 1 T^"' (e) 



'i^)/-T[°\e)T^^\e)/T[^\e)T^''\e) j T[^\e 



(4.6) 



Replacing the T (e) by their asymptotic expansions and expanding the quotient 
leads naturally to transseries p.ip . Hence for linear second-order boundary value 
systems of type (14. ip the ladder transseries template p. II) is the general form. 

5. Nonlinear second-order 2-point boundary value problems 

Given the ladder template for linear second-order 2-point BVP, it is natural to ask 
if a generalisation exists for analogous nonlinear boundary value systems. 

Such systems have been studied asymptotically for decades and form the raison- 
d'etre for many matched asymptotic approaches to boundary layers. Significant 
results were provided by Wasow (1956), O'Malley (1969), Howes (1978). Wasow 
(1956) established conditions under which systems of the form 

eu"{x)=Fi{u,x)u'{x)+F2{u,x), w(0) = a, m(1) = /3, < e « 1, (5.1) 

have absolutely convergent perturbative series on the interval < a; < 1. O'Malley 
(1969) considered a systems approach to the problem, and Howes (1978) used the 
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stability of the e = problem to study the existence of boundary, shock and corner 
layer solutions. 

We shall again restrict ourselves to a situation where only a single boundary 
layer exists, so follow Howes (1978) p. 79 and consider his example (E3): 

eu"{x) + u'{x)u{x) - u{x) = 0, u(0) = a, u{l) = /3, (5.2) 

where 1 < /3 < a + 1. Under these conditions (Howes 1978, p. 83) a boundary layer 
exists near to x = 0. Generalisations of what follows could include, for example, 
larger regions of the (a, (3) plane or e-dependent boundary data. 

The conventional MAE approach to (15.21) generates the composite expansion 

/ (/3_i)tanh(M£)+a\ 

"''-'H""'" ^-'^^'»°-(4^) )'"°"' '"' 

Note that while Mc(0) = a, we again have Uc(l) = /? + O (e^^^^^)/") ^ (i. Clearly, 
when /3 > 1, this error may diminish as e — >■ 0"*". However for larger values of e, or 
small positive /3 — 1 this error may become numerically significant. 
The WKB approximation to the solution takes the form: 



UWKB 



{x) ^([^ + (3-1)+ (i±^_|)(^g-.(/3-i+./2)/A (1^0 (^)) (54) 



Clearly we have wwkb(O) — a, but uwkb(1) — P + O (e"^*^"^/^^/^). Hence, as with 
MAE, the WKB template also fails to satisfy the boundary data at x = 1. 

The WKB approach is just the first two terms of a transseries expansion of the 
general solution of nonlinear ODE which takes the form 

oo oo 

w(x)^^C"7/„(x,e), M„(a;,e)~e-"^(^)/^^a(")e^ ^(0) = 0, (5.5) 

with the constant C determined by the boundary conditions. It will here be ab- 
sorbed into the leading order coefficient of ui, namely Oq . The necessity to in- 
clude integer powers of e'^^^'''^ is caused by the nonlinear terms in the differential 
equation which generate successively higher order exponential scales that must be 
balanced (at each a; 7^ 0), at least in a formal solution. 

However, we see immediately that this template still suffers from the same prob- 
lem as the previous approaches. Specifically the boundary conditions then require; 

5]a(")(0) = 5.oa, a(°)(l) + ^e-"^W/^4")(l)" = "(5,0/3, (5.6) 

n=0 n— 1 

at each r — 0,1,2, For the time being, we put aside the issue of being able to 

satisfy the boundary conditions with an infinite sums of terms. The point is that, 
for F{\) ^ 0, the boundary conditions at x = 1 still leave terms unbalanced at 
O (^e~"-^(^)/'^) for any integer value of n > 0. 
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Inspired by the linear example above, we modify the transseries template as 

CXD OO 

u{x)^J2T.^nH^,e), (5.7) 



uo{x) 



,(0,0) 



r=0 



(x)e^ u(f)(a;,e)^e-("^(^)+P^«)/^ 



,(",p) 



^ — ^min('",p) 



{x)e\ 



(5.8) 

The exponent F{x) again slides between F{0) = and -F'(l) as x runs between 
and 1. The initial assumption of the minimum algebraic order in e, is ry^inin^p) = 0. 
However, we shall see later that even in this simple example, this is not always so. 
In contrast to the single exponent e~^^^'''^ in the linear case, we now have an 
infinite number of exponential scales e""^^^^/*^, n = 0, 1, 2 .... At a; = these all 
degenerate to e""^^"^/*^ — 1 and so can be made to balance at each 0{€^) with the 
corresponding terms from uq. At a; = 1 these scales become e""^^^-'/'^, n = 1, 2 . . . . 
To satisfy the boundary condition at each exponential order thus again requires the 
inclusion of additional series prefactored with scales e"^^^^-'/'^, p = 1,2.... However, 
substitution shows that these boundary series will multiply series in e^''^^^''*^ in 
the uu'{x) term to produce further contributions with prefactors e'-'^^'^^i^)+P^W) / ^ ^ 
ri,p = 0, 1, 2 . . . . All of these series have to be balanced at a; = 0, 1 to satisfy the 
boundary data. Fortunately this is possible, as explained in the diagrams below. 

-(nF(x)+pF{l))/c 






X = 



0<x < 1 



Figure 3. Array structure of exponentials within the double transseries (|5.7j|5.8p . Dots 
joined by lines denote exponential scales e"'"^'^'^''^'^"' "^ that balance at the value of x 
shown. Arrows denote the presence of additional dots corresponding to additional balanc- 
ing exponential scales outside the bounding box of the diagrams. 



From (|5.8p the exponents labelled by integers (n,p) generate not a ladder of 
exponential scales, but a planar lattice (see figure [3]). The heights of the lattice 
points are at the exponential scales e~'^'^^'^^)+P^i'^)) / '^ . At x = 0, 7iF(0) = 0, and so 
balances occur between n = 0, 1, 2, 3 . . . series prefactored by e~P^^'^'^/'^ at each value 
of p. As X increases to 1, F{x) varies and the planar lattice tilts. When cc = 1, the 
prefactors become e~'''^+p)^W/f- and the scales have reconfigured so that balancing 
now takes place between series prefactored by exponentials with n + p — K for 
fixed K (along diagonal rows on the right in figure 3). We shall demonstrate this 
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reconfiguration in detail by explicit calculation for the system (j5.2p in the following 
sections. 

6. Derviation of the coefficients 

Substitution of dH^]) into ([53) and balancing at O (e-("^(^)+P^(i))/<^) gives 

n p 

eui%ix) + E E"™ (^)"-™U^) - ^n\^) - 0. (6.1) 

ni—O q—0 

For fixed {n,p) we can substitute (|5.8I) into (16. ip and balance at 0{e^). We obtain 
the following balances at {n,p) = (0,0), 0{£^) and for n ^ 0, 0{l/e), respectively 

4°-t.(^) + E 4°'°n^)4°:l(^) - 4'^'\x) = 0. (6.2) 

s=0 
n—1 p 

n^F'{x)ai"^P\x) - ^ j](n - m)a[,"'^\a;)4"-"'^-^^(x) = 0. (6.3) 

m=0 s=0 

These recurrence relations now generate the leading order asymptotic terms, 
with associated arbitrary constants Cr'^ . From {n,p) — (0,0) at C(e°) we find: 

ar^(:r)=x + Cr^ (6-4) 

At 0{e) when (n,p) = (0,0), we have 

at'I+-r'at°'+ar'a^^f-ar\^) - ^ a^°'(.) = cf ^°). (6.5) 

It is then simple to prove from {n,p) — (0, 0), 0{e'^), that 

a(o^o)(a;) = C("'°), r>0. (6.6) 

F(a;) is deduced from (n,p) = (1,0) at 0(l/e): 

F'(x)4''*"(:r)-4°'"^(.T)4'-°^(x)=0 ^ F'{x)^a^^'"\x). (6.7) 

This assumes that a^ (x) ^ 0. This result allows us to recast (|6.3I) as a recurrence 
relation to generate Oq ''' (a;). For nonzero values of n at 0{l/e) we have: 

n— 1 p 

n(n-l)ar'4"'^^(-) = ^ E(" " -)«o""^(^)«^'"'''"^^(-) 

m— s — 1 



(m,0) (n-m,p) 
?n— 1 






It is straightforward to show from (j6.8p that 

(n,0) „ (0,0)/ (1,0) ,„ (0,0)\" ,r.n\ 

"o = 2a^ (^a^ V2a^ ^J . (6.9) 
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The term Oq ' (x) can be found, from {n,p) — (1,0), 0{e^) and using (|6.6|) : 

- 4^'")(.) (i + ar\^)cr') - «r"^(-)4!;°^ = o, (e.io) 

^ a^°)(:.) = c(i^°)e-<"'(-+<"'")/(^ + C^'^'°^ (6.11) 

When n = 0, p ^ 0, further terms can be deduced to be constants from the 
0{e^) equations by first observing that, when p — 1: 

(6.12) 
Substitution of this result into the corresponding equation for p > 1 gives: 

f]4°-)(x)4°r^'(-) - «^'^(-) = 0. (6.13) 

s=0 

It follows by iteration that: 

4"'^)(:.)=C^^\ p=l,2,3,.... (6.14) 

We find these constants in the next section. 

7. Boundary conditions reconfiguring the exponential scales 

We now apply the boundary conditions, u(0) — a, u(l) — (3, balanced first at each 
exponential order, and then at each algebraic order of e. The balancing of terms 
within the relevant exponential scales is illustrated graphically in figure 21 

In what follows, we set ar '^ (0) = Cr'^. 

At a; = 0, since F{Q) = the exponential scales O (e~"-^(^y'^^ all collapse to the 
same degenerate level and the distinct exponential scales are given by O (e"^^'-^-'/'^) . 
So for fixed p, the terms Or '*' balance at order O (^e~P^^^'>/'^e^) according to: 

^a("''')(0) = SroSpoa, ^ E«o"'°^(0) = «, ^^^"'''^(O) = 0, {p,r ^ 0). (7.1) 

n=0 n— n—0 

At x = 1, the template (|5.7ll5.8p reconfigures the exponential scales. The orders 
O (e^("+'')^(^)/^) with n+ p — K, K — const, become degenerate. Within these 

scales the terms Ur'^, consequently align at each algebraic 0{e^) according to: 

K K 

^a("^^-")(l) = Sr^SK^P ^ 4"^o)(l) = 5,0/3, ^a("^^-")(l) = 0, (n,r ^ 0). 

n=0 n=0 

(7.2) 
We make progress by first considering the conditions at a; — 1. Using (j6.4p . 
(16. lip . (|7.2p we can identify immediately that 

ar^(l) = l + Cr^^/3 ^ a'^r\x)=x + p-l, (7.3) 

a(o.o)(l) = 5,0/3 = C^")=0 r>0, ^ ^''"^(x) = C^''°V(^ + /? " 1)- (7-4) 
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4''''(1) 
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X = a; = 1 

Figure 4. Balancing of terms ar'^ at fixed r between transseries with varying (n,p) 
(|5.7j|5.8|l at a; = and x — 1 due to the realignment of the exponential scales. Terms 
joined by the grey lines balance at the values of x shown. At a: = at each fixed p, an 
infinite sum over all terms n must satisfy the boundary conditions at the left hand end 
of the interval. At a; = 1 terms with n + p = K for each fixed K = 0, 1, 2 . . . must be 
summed to satisfy the boundary conditions at the right-hand end of the interval. 



Hence, returning to x = we have, for example, 

oo 



(7.5) 



Using the form (j6.9p we may smn this infinite series formally as a geometric pro 
gression and obtain 

2a^")(0)a^°)(0) 



,(0:0) 



(1,0). 



2ar'i0)-ar>i0) 



= a-^+l 



a-/3 + l, 



2(/3-l)C<^'")(0) 
2(/3-l)2-C^i^°)(0) 



Returning back to a; = 1 we then have: 



a 



(1,0) ^ 2(a-/3 + l)(/3-l)^ 
{a + (3-1) 



(7.6) 



(7.7) 



al^^'\l)+ai'''\l)^0 ^ Ct'^ + ^ 



(1,0) 



^ C, 



(0,1) _ 2(a~/3 + l)(/3-l)^ 



/3(« + /3-l) 

(7.8) 
Hence we finally have the leading orders of the first three {n,p) transseries as: 



,(0,0) 



(1,0)/ ^ _ 2{a - (3 + 1){I3 - If „(o,i)/ ^ _ „(i,o). 



From (|6.14p and (j7.2p we can also deduce that for this example ar (x) — 0, r > 1. 
Note that this means that for the special case of (|5.2p . the expansion Uq ' will 
not generate a Stokes phenomenon, since it truncates after the first term. This also 
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reflects the fact that there are no internal boundary layers. Note also this means 
that the exponent F{x) is not here directly connected to a factorial-over-power 
late-term ansatz (Dingle 1973; Berry 1989; Chapman et al 1998). 

We may resum the n-transseries for p = using (17. 9|) and obtain the following: 



«tra„s(x) = 4"^°)(x) + Y. e-"^("^/^4"^°)(a:) (1 -I- 0{e)) 



(7.10) 



n=l 



JO.O) 



(x) 



24"-°^(x)a^°)(x)e-^(-)A 
24°'°)(a;)-4''°\a:)e-^(-)/^ 



(l + 0(e)).(7.11) 



For most practical purposes, this result may be all that is needed. 

A comparison of the numerical errors for the MAE (|5.3|) . WKB (15. 4p and resum- 
mation (|7.10p is shown in figure [S] for typical values of e = 1/10 and 1 with fixed 
a = 3/2 and /3 = 2. First, for small e (left hand graph), although the errors com- 
mitted by the approximations at a; = 1 are too small to display on the scale of the 
graph, the transseries approximation is uniformly numerically better than either 
the matched or WKB solutions. Secondly, for larger values of e both the WKB and 
transseries approximations are uniformly better approximations than the MAE. 
The latter violates the boundary condition at a; = 1 by a noticeably larger amount. 
Although it satisfies the boundary condition at x = 0, the transseries approximation 
is not uniformly better than the WKB near there, agreeing more with the MAE. 
However the transseries approximation quickly does beat WKB as x increases and 
is definitely better at a; = 1. Note that the transseries approximation (|7.10p only 
contains terms that have p = 0. If we summed the n-series with p = \ and added 
these in at 0(e), the error at a; = 1 would diminish. In the next section we show how 
the higher order balances reveal unexpected behaviour in the nonlinear template. 







X 



0.2 0.4 0.6 0.8 1.0 

e= 1/10 
Figure 5. Comparison of the percentage relative errors, 100|1 — u. 
matched (dotted), WKB (dashed) 
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8. Higher order balancing 

Moving on to higher exponential orders, at a; = 1 we need to balance terms involving 
the (n,p) = (2,0), (1,1), (0,2) series (cf. figure i]): 



,(0,2) 



(1) 



,(14) 



{l)+af'\l)^Q. 



5.1) 
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The behaviour of the {n,p) — (2,0), (0,2) terms have already been derived above: 



(ap •'(x)j 01'^, _ fl ^ n2^fl _ n4 



(2,0), , _ r" ^^^; ^ 2{a-p + l)\p-lf (0,2), , _ ^(0,2) 

2ar^(x) (.T + /3-l)3(a + /3-l)2 



"" ' ^""^ - ~jm7Zr - TTTJ^WuTTW^T)^' "" ' (2^) - c^o ' ■ (8-2) 



Seeking the aj, ' (a;) term we use (|6.1I) with (n,p) = (1, 1): 

-i!L + 4"'4:i + 4^^4°i + "^4^1 + 4^^<i - u['^ = 0. (8.3) 

Assuming r^in{l, 1) = (cf. (IS-SI) ). at 0{l/e) in (|8.3p . after simphfication we have 

F'(x)4°''^(:e)4''°Hx)=0. (8.4) 

This obviously leads to a contradiction, since prior calculations show that neither 
of the conditions F'{x) = 0, Oq ' '{x) — 0, or Oq '(x) — hold. 

We may resolve this impasse by taking ri„in(l, 1) = —1. For then, at 0{l/e^), 
(ESI) is satisfied by F'{x) = a[,°^"^ and at 0(l/e) we have 

ar\-)a^'i%) + a^-,'\-) = -ar^iM'^'H-h^'i-), (8-5) 

with general solution, 

(1,1), , _ Ciy^/3(a + (3-ir + Ax{(3 - l^a ~ P + 1)' ,_. 

"-' ^'''~ P{x + P-l){a + (i-lY ■ ^^-^^ 

Given that, from above, there is no term a_{ (0), the boundary condition that this 
term must satisify at x = is (cf. figure |3]), 

5^aLr^(0) = 0. (8.7) 

r(=l 

By examination of the ©(l/e^) balances in the (n,p) = (l,p) equations (|6.ip . it is 
possible to show that 

aL"i'^(a;) - n (4^'")(a;)/24"'°)(a;))""' aLV^(^)- (8-8) 

Hence we can sum (|8.7I) to obtain (cf. (18. 6p ) 

4aLV'(0)(a^°'(0)) ^ d Df ^ = 4a;(/3 - l)4(a - /3 + 1)^ 

(2«r(0)-ar(0))^^ ^ ""^ ^'^ /3(. + /3-l)(. + /3-l)- 

(8-9) 
Consequently, since a_{ {x) = we can now determine the constant a_[ = C_{ 
from the boundary condition at a; = 1: 

aLY^(l) +aLY^(l) + aLY^(l)=0 

^ aLY^(.) =-aLV^(l)^-li^^ll%^il±ll!. 



(8.10) 
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Finally, it is then possible to resum the contributions from a_-l to obtain 
- .(".Dm 4a<Y^(a:)fa|,''°^(x)Ve-(^(-)+^(i))/^ 

^0 ^ e(2ar)(x)-a(^''')(x)e-^WA)' 

The exponential asymptotic approach has revealed the unexpected and counterin- 
tuitive presence of 1/e terms in the expansion as e — >■ 0"*". However, these 1/e terms 
are prefactored by e~^(^^/' with F{1) > 0. Hence they are not only beyond all 
algebraic orders, but also do vanish exponentially fast as e — >■ 0+. 

The existence of an 0(l/e) term in the {n,p) = (1, 1) series forces the presence 
of 0{l/eP) terms in series with p > 1. The corrected form of the template (j5.8p has: 

"-''^("'^) = { -Floor[(n+p)/2], n < p. ^^'^^^ 

As above, since the negative powers of e are multiplied by e^^^^^''"^, they will remain 
technically beyond all algebraic orders and vanish as e ^- 0+ . 

The presence of these terms means that in generic cases, as p increases, addi- 
tional work must be done to obtain resummed expansions accurate to 
Q(^Q-n{F{x)+F{i))/e^y For example, it is a simple algebraic exercise to show that the 
calculation to derive the resummation analogous to (|8.1ip of the 0(6°) terms with 
p=l, 

^g-«(F(.)+F(l))/.^(«,l)(^)^0^ (8.13) 

actually requires first a complete derivation of all the terms in a^"' (x). 
A detailed set of calculations generates the following results for n > 1: 

a[--'\x) = n (a^°V2a^°')"~' a['-'^ + Kja^^ (a^°V24"'°')" , (8.14) 

+nin - l)(4i'°))""' / {2ar^y~\['''^a^^i'^ 

-2(n-l)(a^°V2ar)"«^'^ 

+2i„(4^'"^)""/(2a^°^)"^^LV'. (8.15) 

The coefficients i^„ in (j8.14p and L„ in (j8.15p satisfy the recurrence relations 



(n-l)i^„/2 = ^i^„-(2n+l)(n-l)/n, (8.16) 

n—1 n— 1 

(n-l)L„/2 = "^L^+"^{n-m)K,n-{2n+l){n-l). (8.17) 



m—2 771—2 



Article submitted to Royal Society 



n=2 



Exponential asymptotics and boundary value problems 17 

A z-transform of these relations generates 

";.,r.Hm. ^";-j^;;'"' + ^" + ;'^f-" . (8.18) 

f:u."^Gw. ^'"";i%f''"» +g^^. (S.19) 

where Li2(t) is the dilogarithm function (DLMF 2010). 

The sum over n of (|8.14p at cc = can be achieved using (j8.18p and according 
to the boundary value, be set equal to zero. It only depends on unknown a\ ' and 
the already known Oq ' and aj, . Coupled with the general solution of 

(0,0) (1,0) (1,0) _ (1,0) ,^ „^s^ 

^0 ai,x +^1 -H,xx^ (.»./!0j 

we can use this sum to determine a\ ' {x). This can be substituted into (j8.15p . 
which, in turn can be summed at a; = using (|8.19l) to find aj, ' (a;). In turn (I8.15P 
can then be multiplied by e"'"^'^^-'^^'^"'^^^/^ and summed over n to obtain (|8.13p . 
The resulting expression is large and so this is left as an exercise for the reader. 

In principle this scheme can be repeated for each p. For given p, the correspond- 
ing 0{e) estimate will involve a sum over n of Oq '^ (a;)e"'^"^^^^+P^'^^)). To deter- 
mine the Oq '^ (x) will likely require the determination of all terms from a,."''' (a;), 
> r > r-[ain{n,p) beforehand. This will become an increasingly arduous task. 

However, as the transseries template contains three sums (in n, p and r) there 
is the potential to alter the order of resummation of the terms to improve the 
numerical accuracy of the approximation. 

9. Alternative resummations 

The above resummations are over n for each p and fixed order e^ . Alternatively it 
may be possible to reorder the resummations, for example resumming first over p 
and obtain numerically better agreement with the exact solution. 

Due to the particular form of (J5.2I) we may proceed here down this route by 
reverting to a single transseries formulation of the form 

oo oo 

u{x)=Y,\''u^{x-t), So(a;,e)=^a(°)(a:)e^ (9.1) 

fe=0 r=0 

where A is here just an ordering parameter that will be eventually set equal to one. 
In the absence of the boundary data, the Uk{x;e) can be regarded as providing 
terms of order O (exp(— fcF(x)/e)). Substitution into (|5.2p and balancing at order 
A'' generates at ©(A"): 

tui^{x) + uq{x)uq(x) — uo{x) = (9-2) 

and balancing at 0{e^) and setting a}. (1) — 6rol3 we have 

dl^\x) = Sroix + 13 - 1) uo{x)=x + ^~l, (9.3) 
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as above. Note that the truncated nature of uo{x) is due to the particular system 
under study. In general it would be only expressible as a formal, divergent series. 
At order 0{X) we then have 



eu'({x) + {x + j3 — l)ui{x) — 
Applying the boundary condition, e.g., that (with A = 1) 



0. (9.4) 



uo(0)+Mi(0)=a, wn(l)+wi(l)=/3, 



(9.5) 



we find 



ui{x) = {a- 13+1) erf 



x + 13- 1 



— erf 



/erf 



/3-1 



erf 



(9.6) 
The numerical approximation of uq{x) + ui{x) is better than any of the MAE or 
resummations so far (see figures [5l ^. This is because the boundary conditions 
can here be satisfied at both ends by the leading two orders uo{x), ■ui{x). Note 
that the exact forms (j9.3p . (j9.6p must incorporate resummations of the r series 
in algebraic powers of e, and this is here actually the reason why both boundary 
conditions have here been able to be satisfied. That we can achieve this is due to 
the fortunate peculiar form of (|9.4[) due to the truncation of the expansion in (|9.3p . 
This serendipity would not extend to general systems and the r sum would not 
have been automatically resummed. 

We can also solve the linear inhomogeneous equation at C'(A^) for (|5.2p with 
boundary conditions M2(0) = ^2(1) = 0. Addition of M2(a;) to uo{x)+ui{x) improves 
the numerical agreement even further (figure [5]). 




e= 1/10 



e= 1 



Figure 6. Comparison of relative percentage errors for e = 1/10, 1 with a — 3/2, /3 = 2 
for resummed approximations using only uo{x), uo{x) + ui{x) and uo(x) + ui{x) + U2{x). 
Note the satisfaction of the boundary condition at both ends and the rapid apparent 
convergence. 



The form of ui{x) in relation to the double-transseries template requires an 
explanation. If we substitute the asymptotic expansion erf(s) ~ 1— exp (— s^) /\/t^s, 
s — > +00, into ui{x) we find 



ui{x) ^ v{x) — v{l) v{x) = 



(/3-l)(a-/3 + l)e--P'(^)/<^ 
(:e + /3 - 1) (1 - (/? - l)e-^(i)/V/3) ' 



(9.7) 
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where F{x) is as defined in the sections above. At first sight, the form of v{x) 
suggests that ui{x) might be interpreted as the ©(e") term of the resummed p- 
series for n = 1, (|7.10p . However a simple calculation shows this not to be the case. 
In fact ui{x) does not satisfy the same boundary condition as the resummed (j7.10p . 
The reordering is more complicated with terms of 0{e~^^'''^) appearing also in the 
higher order u^ ■ Higher order itu satisfy inhomogeneous hnear differential equations 

eu'l{x) + {x + (i-l)u'f.{x) = /(uo,"o'"i''"i>---"fc-i)' "fc(0) ^ ""^(1) = 0. (9.8) 

The solutions of the homogeneous form equation will be of the same form for all k: 

OQ 

A^{t)v^{x)+Bt,{e)v2{x)e-P^-y\ «,(:e) ^ ^ 4')(x)e^ / = 1, 2. (9.9) 

Terms of C'(exp(— ji<"(a;)/e)), < j < k will arise from the inhomogeneous terms 
and so the Ak{e) and Bk{e) will contain terms in 0(e~-'^'^^^/'), < j < fc to 
satisfy the boundary conditions. Hence the expansion (j9.ip can be seen as non- 
trivial reordered sum over p and r, ordered by the index n. For this particular 
system, due to the truncated nature of mq, we are able to solve for the Uk exactly 
and so effectively achieve this double resummation. For situations where uq does 
not truncate this will not be achievable in general. 

Note that the transseries template is forced by the presence of boundary values, 
but the form of the lambda series is driven by the differential equations. It is an 
open question as to what is the ordering of terms in the transseries that optimises 
the numerical agreement with the exact solution, when only a finite number of 
terms is taken in the approximation. 

10. Discussion 

We conclude with a brief comparison of our nonlinear results with those of Wasow 
(1956) and a discussion of further work. In our notation, Wasow considered (|5.ip 
with Fi{x^u,e) and F2(x,u,e) regular analytic with respect to u, e and C*-^^ in x 
in a region of (x,u, e) space containing u — uq(x), < a; < 1 and e = 0. Wasow's 
expansion takes the form 

oo 

u{x,e) — 2_\vr{x,e)^^ , Vr{x,e) = e~^'^^'''^Wr{x,e), r > 0, 

where fx — a — fo(0, e) — a — /3 + I, v{l, e) = /?, the Wr{x, e) are bounded and F(x) 
is as defined above. He proved that the u,. expansion is uniformly and absolutely 
convergent with respect to e, /i and x for 0<e<ei, |a — /? + l|<Mi,0<a;<l. 
Wasow then reordered the terms to give 

oo 

u{x,e)^ ^c„(a:,a,e)e-™^(^)/^ 

m=0 

as uniformly and absolutely convergent. He stated that the coefficients Cmix,a,e) 
are regular analytic and possess a power series in e. He omited the full proof as it 
was, in his own words, "somewhat detailed". 



Article submitted to Royal Society 



20 C.J. Howls 



Wasow's results are existence proofs. He took the resumraed route directly, his 
"reordering" is related to the transseries. Here we have taken the reverse route of 
solving for the transseries first and then resumming as we have been motivated by an 
exponential asymptotics approach. Our explicit calculation reveals the intricate and 
subtle nature of the form of Wasow's Cm or Wr coefficients and reveals an interesting 
realignment of exponential scales that are required to satisfy the boundary values. 
The template has been dictated by the boundary values and is independent of 
the equation and so is expected to be more generally valid. Wasow's statements 
that the coefficients "possess an asymptotic expansion in powers of e" is strictly 
incomplete, since it ignores the presence of series beyond all orders, and in inverse 
powers of the asymptotic parameter. Given all this, Wasow's work suggests that 
our resummations may be convergent, at least for small values oi jj, = a — /3 + 1. 

The preliminary work of this paper has opened up a multitude of potential 
follow-up problems in exponential asymptotics. Further work could include: rigor- 
ous justification of the exponential approach, with a proof (or otherwise) of conver- 
gence; derivation of proper error bounds for the reordered summations; extensions 
to higher order BVPs; adjustment of the template to include interior layers and 
shocks; a full hyperasymptotic treatment based on the fundamental transseries 
template; extensions of the transeries template to PDE BVPs. Some of these will 
be discussed elsewhere. 

The author acknowledges the hospitality of the Pacific Institute for the Mathematical 
Sciences and the Department of Mathematics at UBC, Vancouver, where some of this 
work was undertaken, and AB Olde Daalhuis for very helpful discussions. 

Appendix A. Multiple-scales exponential asymptotics 



Motivated by the form of (j3.7p we seek, a priori, a multiple-scales solution of (12. ip 
that corresponds to the summation of all the p-exponentials in (13.11) . We use the 
ansatz: 

oo oo 

u{x; e) ^ ^ Wr{x, Xy + e'^^^^/' ^ Vn{x, V)e^ (A 1) 

where the scaled "variable" is the constant X = F{l)/e. 

We substitute dSH) into ^^ and balance at O^) and ©(e^^^^^/'^e''), ignoring 
the e-dependence in the X terms. Since X is a constant, it does not actually generate 
a derivative in d/dX. We thus obtain recurrence relations that are identical to those 
we obtained when we substituted (j2.5|) into (12. ip . Hence the recurrence relations 
for W and V are, for r > 0, 

Wr{x,X) = (c,-W;_i(x,V))/(2x + l), V;{x,X) = {V;'_i{x,X))/{2x + l), (A2) 

with W^i{x,X) — V-i{x,X) = and prime denoting x-differentiation. 

The constants Cr are determined from a modified set of boundary conditions. 
The terms in the expansion in (J2.10I) , (|2.1ip and (j3.4p have been balanced at orders 
of 0(e~^'^-'/'^e''). Starting from (jA ip due to the apparent similarity of (|2.5p and 
the second row of (jA II) . we find that we go around in a circle. The second sum in 
(lA ip cannot satisfy the boundary condition at a; = 1 exactly, without including 
further series, for the reasons outlined above. Thus the boundary conditions have 
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to be modified. Instead of (l2.10ll2.lTt we liave: 

Wr{Q,X)+Vr{Q,X)=5rOa, Wr{l,X)+C-^Vr{l,X) = 5rQl3. (A3) 

Note tfiat the exponentiai in (jA3l) is actuaily e"^^"'^^/'^, precisely tlie order of ne- 
glected terms in (|2.1ip which lead to the ansatz p. II) . That we can include this 
exponential term now is because we treat X as varying on a different scale to e. 
From (jASp the Cr in (IA2p can now be found. A short calculation gives 



(T4^;_i(l, X) - WUJx, X)) - 3c-^(T4^;_i(0, X) - WUJx, X)) 

^^("'^) = (l-3e-^)(2.x + l) ' 

Vr{x,X) = (M^;_i(0)-M^;_i(0))/(l-3e-^). (A 4) 

with initial terms Wo (x,X) = 3(/3- ae"^)/ {(1 - 3e"^)(2x + 1)} imAVo{x,X) = 
{a — 3/3)/(l — Se"'^). Inserting the value of eX = F{\) = 2, to leading order in e 
we recover (I3.7p . Note that by neglecting the terms in e""^ these relations reduce 
to those of (j2. 12112.1 5)1 ■ as they should. 
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